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The process for CO 2 removal from flue gases was modeled with RateFrac. It consists 
of an absorber, a stripper, and a cross heat exchanger. The solvent used in the model 
contains about 30 wt % monoethanolamine (MEA) in water. MEA reacts with CO2 in 
the packed absorber The finite reaction rate requires a kinetic characterization. The 
RateFrac absorber model was integrated with a FORTRAN user kinetic subroutine to 
make the model consistent with the interface pseudo-first-order model and with a re- 
gressed Electrolyte-NRTL equilibrium model. It was adjusted with laboratory wetted 
wall column data and field data from a commercial plant. Sensitivity analyses were 
performed on process variables to find operating conditions at low steam requirement. 
Many variables strongly affect the process performance, but an overall optimization 
shows that there are no economical ways to reduce the steam requirements by more 
than 10%. The reboiler duty can be reduced from that of a base case representing 
current industrial operating conditions, by 5% if acids are added to the solvent, by 10% 
if the absorber height is increased by 20%, and by 4% if the absorber is intercooled 
with a duty of one-third of the reboiler duty. The power plant lost work is affected by 
varying stripper pressure, but not significantly, so any convenient pressure can be chosen 
to operate the stripper. 



Introduction 

The removal and sequestration of COj from combustion 
gases is an important technological alternative to address 
global climate change. Absorption/stripping with aqueous 
monoethanolamine (MEA) has been commercially applied in 
small plants for CO 2 recovery. Aqueous MEA is an effective 
solvent for CO 2 capture, but a system for 90% CO 2 removal 
can reduce the efficiency of a power plant from 40% to 30% 
(lEA Greenhouse R&D Programme, 2002). 

The absorption/stripping process is shown in Figure 1. Flue 
gas contacts the aqueous solvent at 1 atm and 40-60°C, in a 
countercurrent, packed absorber. Typical sources of flue gas 
include gas-fired turbines, giving 3 mol % CO, and coal-fired 
plants, giving 10-12% CO.. 

One version of this process uses approximately 30 wt % 
MEA; another uses 15-20 wt % MEA (Liljedahl et al., 2001). 
Typically the lean solvent has a CO2 loading of 0.1-0.2 
mol/mol MEA, and the rich solvent has a loading of 0.4-0.5. 



Correspondence concerning this article should be addressed to G. Rochelle. 



The rich solvent is stripped by steam in a countercurrent, 
reboiled column at 1.5-2 atm and 100-120°C to produce pure 
COj. Heat is recovered from the hot lean solvent by cross- 
exchange with cold rich solvent. The lean solvent is typically 
cooled further to approximately 40°C. 

Integrated models for this process have been created; they 
either use commercial software or language codes. TSWEET 
has been used since the early 1980s for modeling of acid gas 
removal (Holmes et al, 1984). Another software package in 
use is AMSIM, which uses a rigorous nonequilibrium-stage 
model (Zhang et al., 1996). Weiland and Dingman (2001) de- 
veloped a program called ProTreat for the rate -based simula- 
tion of columns. Along with commercial packages, programs 
written in Fortran or Visual Basic have been developed. These 
programs have the advantage of being specific for amine gas 
treating. A model by Al-Baghli et al. (2001) uses this method. 
This approach is usually slower, and presents challenges in 
simulating the whole process. 

The purpose of this work is to understand how the design 
variables affect each other at the level of the whole process. 
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In addition to the species included in the Austgen frame- 
work, this model also included heat stable salts (HSS) mod- 
eled as formates, formed according to the irreversible reac- 



Figure 1. Absorption /stripping process for CO2 re- 
moval with allonolamines. 



not of the single absorber and stripper. Aspen Plus was found 
to be suitable for this type of analysis, and it was chosen as a 
framework for the model built in this work. The model was 
built to analyze the effect on energy requirements of several 
process variables. The goal was to find operating conditions 
that allow CO 2 removal with less energy. 

Model Development 

The absorber and stripper were both modeled with Rate- 
Frac, a rate-based model framework in Aspen Plus. In the 
absorber, the reactions involving COj were described with a 
kinetic model. In the stripper all the reactions were set to 
equilibrium, due to the higher operating temperature. This 
approach required that a rigorous thermodynamic model and 
a rigorous rate model be implemented in Aspen Plus. 

Thermodynamic model 

The thermodynamics were described with the Electrolyte- 
NRTL framework, developed by Chen et al. (1979), and mod- 
ified by Mock et al. (1986) for mixed-solvent systems. The 

model reproduces the one developed by Austgen (1989), with 
some interaction parameters regressed in order to match the 
data of Jou et al. (1995), considered more accurate than those 
of Lee et al. (1976) and others, used by Austgen. The reac- 
tions incorporated in the model are the same as those in- 
cluded by Austgen (1989). The values of the constants and 
the interaction parameters can be found in Freguia (2002) 



MEACOO" -l-HjO^-^MEA-l-HCOj (la) 

C02-I-2H20«^^HC03" -FH30+ (lb) 

HCOa' +H20 <^°"''°''> C03~~ +H3O+ (Ic) 

MEAH + + H2O <^°™^"> MEA+ H3O + (Id) 

2H20«^H30 + -)-OH- (le) 



MEA + HCOOH ^ MEAH + HCOO - 



(2) 



HSS are normally present in absorption/stripping processes 
for acid gas removal. They are products of MEA neutraliza- 
tion with strong acids. Typical forms of HSS are sulfates or 
formates. The HSS were modeled in Aspen Plus by including 
formate as a component and assigning to it the same interac- 
tion parameters as MEACOO". No VLE data were available 
with formate in the system. The important properly of the 
HCOO" anion is its charge. This extra charge has to be bal- 
anced by MEAH"^, thus making some of the MEA less avail- 
able for CO2 absorption. The presence of HSS alters the 
equilibrium, increasing Pqq^ at a given loading, which is re- 
defined according to Eq. 3. This fact causes the absorber 
driving forces to be reduced, and the driving forces in the 
stripper to be increased 

[CO2L, 



... — ^ - " - [MEA].„.-[HSS] - [MEA];i 
The concentration of HSS is also given by the loading 
[HSS] 



Rate model 

The rate model used in the absorber was made consistent 
with the Electrolyte -NRTL thermodynamic model and with 
the interface-pseudo-first-order (IPFO) approximate model 
for mass transfer, through a Fortran kinetic subroutine for 
RateFrac. In the rate model the reactions were the same as 
in the equilibrium model, except reactions la and lb, which 
were replaced by kinetic reversible reactions 5a and 5b. The 
direct reactions were given second-order rate expressions, 
given by Eqs. 6 and 7 

C02 + MEA-hH;0-^2^MEACOO- -hH30+ (5a) 



-"COj - MEA - ''2,MEA*-t0f*MEA"CO2 
RcOo - OH - = ^2 OH - CLXq^ - flco. 



(6) 
(7) 



Reaction 5b is almost always negligible in the reaction 
boundaiy layer, but it is necessary to provide a route for the 
reversible production of bicarbonate. 

The IPFO model assumes that a consequence of the fast 
reaction kinetics of COj is that the concentration of eveiy 
other species in solution is constant at its interface value. This 
assumption has the advantage that the COj flux can be cal- 
culated analytically, using Eq. 8, where the driving force is 
expressed in terms of activities, calculated using the Elec- 
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trolyte-NRTL model 

Nco, 




( ^2,MEA[MEA]i + yfcj OH - [OH - ],) — ^ 



X(«co„,-«So„,) (8) 

The term flgg^ is defined by Eq. 9 

^2,MEA '^MEACOO^ ^H30+ ^ ^2,OH~ ^HCOa' 
^ -^^COz-MEA '^HjOTmEA -^^COj-OH- ')'oH " 

''^2,MEA-*MEA ^ ''Z.OH" -^OH" 

(9) 

In the equation the rates of reverse reactions 5a and 5b 
were calculated using their equilibrium constants. Equation 9 
is an arbitrary allocation of the production rate of COj be- 
tween reactions 5a and 5b. 

Equation 10 defines the rate constant used for reaction 5a. 
The temperature dependence is derived from wetted wall 
column data from Dang (2001). The constant A was adjusted 
to match commercial plant data, provided by Fluor Daniel, 
Inc. Equation 10 with Dang's data implies rate constants 
higher than the values measured by Hikita et al. (1977) by 
factors of 2 to 5 at 40-60°C. This difference is significant, 
and it can be attributed to the increased solution ionic 
strength in loaded solutions. This is consistent with the re- 
sults of Pinsent et al. (1956), who showed that ionic strength 
as 1 M NaCl increased the reaction rate of COj with ammo- 
nia by a factor of 3 

904.6 

lnfc2,MEA=^ — (10) 

Equation 11 gives the expression for the rate constant of 
reaction 5b, from Sherwood et al. (1975) 

6658 

In^oH- =31-396 — (11) 

The diffusivity of COj was calculated using the Stokes- 
Einstein relationship proposed by Pacheco (1998), given by 
Eq. 12 

(^) '''' 



The viscosity of the solution and that of water were calcu- 
lated using the expression developed by Weiland (1996), 
which accounts for temperature, MEA concentration, and 
CO 2 loading. The addition of MEA to water increases the 
viscosity, thus reducing the diffusivity. 

The activity coefficient of COj (7002) serves to transform 
the concentration-based diffusivity into an activity-based dif- 
fusivity. It is calculated by the Electrolyte-NRTL model. 

Other RateFrac methods 

The RateFrac module of Aspen Plus was used to integrate 
the heat transfer and multicomponent mass-transfer relation- 
ships for a packed absorber and stripper. Among the many 
rate phenomena included, the most important were COj liq- 
uid film mass transfer, HjO gas film mass transfer, and gas 
film heat transfer. This model used the rigorous Electrolyte- 
NRTL thermodynamics to predict enthalpies, equilibrium va- 
por pressures of COj and HjO, and solution speciation. 
Gas-film and liquid-film mass- and heat-transfer coefficients 
were obtained as a function of liquid and gas flow rates and 
properties with the default correlations of RateFrac. 

Model validation 

The rate and equilibrium model were used to simulate a 
base case for the absorption/stripping process. The base case 
was developed on proprietary data from a commercial plant. 
CO2 loading, packing height, and stream flow rates will not 
be disclosed here. The reboiler duties are normalized to the 
moles of CO2 removed and they are divided by a typical op- 
erating value. 

The constant A in Eq. 10 was adjusted in order to match 
the base-case removal. Table 1 compares the plant data with 
the model results, for the base case and three other test runs. 

In the base-case simulation COj removal is well matched; 
the reboiler duty is overpredicted by 3.5%. In the other three 
runs the main mismatch is in the reboiler duty, which is over- 
predicted by 10 to 30%; the COj removal is underpredicted 
by up to 4%. 

Effect of Process Variables on Energy 
Requirements 

Model methods and base-case conditions 

The following conditions were used in the base case, and 
were kept at those values in all the runs, unless otherwise 
stated. 

• MEA concentration, [MEA],,,,: 33.5 wt % 

• HSS loading: 0.1 mol/mol MEA,„, 

• Stripper bottom pressure: 1.7 atm 

• Absorber-solvent inlet temperature: 40°C 



Table 1. Summary of Results of Model Validation with Commercial Plant Data 





Base 
Case 


Model 


Testl 


Model 


Test 2 


Model 


Test 3 


Model 


Feed CO2 (%) 


3.13 




2.87 




3 




2.86 




Vent CO2 (%) 


0.49 


0.49 


0.34 


0.342 


0.56 


0.63 


0.3 


0.38 


CO2 recovery (%) 


85.0 


85.0 


88.0 


87.9 


81.5 


78.6 


89.7 


85.5 


MEA (wt %) 


30 




28.2 




26.6 




29 




Relative 


0.99 


1.035 


1.125 


1.515 


1.02 


1.2 


0.995 


1.135 
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Table 2a. Detailed Simulation Results 



Table 2d. Detailed Simulation Results (continued) 



Solvent Rate: 3%C02, 85% Removal, 33.5 w 
MEA,„„ HSS Loading = 0.1, PfeboUer = 1-77 atm, 

■Zbase case = U -Zstripper/Zbase case = 1 





Lean 


Rich 






Loading 


Loading 


Crel 


0.633 


0.11 


0.417 


1.490 


0.684 


0.13 


0.415 


1.269 


0.735 


0.15 


0.413 


1.145 


0.780 


0.16 


0.411 


1.085 


0.863 


0.18 


0.408 


1.050 


0.912 


0.19 


0.406 


1.053 


0.968 


0.20 


0.404 


1.063 


1.032 


0.21 


0.402 


1.078 


1.299 


0.24 


0.393 


1.143 


1.775 


0.27 


0.383 


1.255 



Table 2b. Detailed Simulation Results (continued) 



Absorber H: 33.5wt. %MEA„„Z,,,ipp,,/Zb.,,e case = 1, 
HSS Loading = 0.1, P.^boiier = 1-77 atm. Optimum L/G 



CO2 in 
Flue Gas 
(mol %) 


C02 

Removal 
(%) 


■Zbase cas 


L/Gmass 


Loading 


Rich 
Loading 




3 


85 


0.6 


3.29 


0.26 


0.321 


1.87 


3 


85 


1.6 


0.82 


0.20 


0.440 


0.95 


3 


90 


0.8 


3.37 


0.28 


0.343 


1.69 


3 


90 


1.0 


1.18 


0.19 


0.368 


1.25 


3 


90 


1.6 


0.85 


0.19 


0.435 


0.97 


10 


90 


0.8 


5.42 


0.29 


0.417 


1.11 


10 


90 


1.0 


2.68 


0.19 


0.442 


0.96 


10 


90 


1.6 


2.61 


0.20 


0.460 


0.90 



Table 2c. Detailed Simulation Results (continued) 

Stripper H: 3%C02, 85% Removal, 33.5 wt. % 
MEA,o„Z3bsorber2/basecase = L HSS Loading = 0.1, 
■Preboiier = 1-77 atm. Optimum L/G 



^strippsr 






Rich 














0.285 


1.174 


0.23 


0.399 


1.142 


0.428 


1.031 


0.21 


0.402 


1.110 


0.714 


0.967 


0.20 


0.404 


1.079 



• Flue-gas water-saturation temperature: 115°C 

• Flue -gas absorber-inlet temperature: 63°C 

• Cross-exchanger hot-end temperature approach to equi- 
librium: 11°C 

The cross heat exchanger was modeled in Aspen Plus with 
two separate HEATERS, connected by CALCULATOR 
blocks in such a way that the temperature approach ol the 
hot end is fixed at ITC and the heat duties ol the two heaters 
are matched. 

In the absorber the packing is modeled with 20 RateFrac 
segments; in the stripper 19 segments represent the packing, 
the bottom segment being the equilibrium reboiler. Some cal- 
culations were performed with 30 RateFrac segments repre- 
senting the same total height of packing. There was less than 
1% difference in the calculated reboiler duty. The top seg- 
ment models a water wash. The rich feed enters the column 
"above" the second segment from the top. A condenser is 



Heat Stable Salts, MEA,„,: 3%C02,85% Removal, 

■Zabsorber/^base case = U Z,,,^^JZ^,^^ = 1, P^eboiler 

= 1.77 atm. Optimum L/G 



MEAtot 
(wt%) 


HSS 
Loading 


V^mass 


Lean 
Loading 


Rich 
Loading 


Grel 


20 


0 


1.25 


0.22 


0.451 


1.201 


20 


0.05 


1.279 


0.21 


0.448 


1.197 


20 


0.10 


1.321 


0.20 


0.444 


1.196 


20 


0.15 


1.314 


0.18 


0.441 


1.195 


20 


0.20 


1.319 


0.16 


0.437 


1.198 


20 


0.25 


1.39 


0.15 


0.432 


1.203 


30 


0 


0.97 


0.23 


0.431 


1.093 


30 


0.05 


0.951 


0.21 


0.427 


1.084 


30 


0.10 


0.99 


0.2 


0.421 


1.074 


30 


0.15 


0.988 


0.18 


0.416 


1.067 


30 


0.20 


0.996 


0.16 


0.41 


1.063 


30 


0.25 


1.021 


0.14 


0.402 


1.066 


36 


0 


0.928 


0.23 


0.407 


1.087 


36 


0.05 


0.86 


0.2 


0.402 


1.069 


36 


0.10 


0.896 


0.19 


0.396 


1.053 


36 


0.15 


0.951 


0.18 


0.387 


1.043 


36 


0.20 


0.871 


0.14 


0.381 


1.04 


36 


0.25 


0.894 


0.12 


0.372 


1.033 



Table 2e. Detailed Simulation Results ( continued) 
Stripper P: 3% CO2, 85% Removal, 33.5 wt. % 

MEA,„„ HSS Loading = 0.1, Z.bsorber/Zbase case = 1> -Z.^ipp,,/ 

Zbasecase = 1' Optimum L/G 

Lean Rich Relative 
^rebouer ^-A^mass Loading Loading j2rei Lost Work 



1.529 
0.824 
0.824 
0.765 
0.765 
0.765 



0.28 
0.17 
0.16 
0.14 
0.13 
0.12 



0.241 
0.243 
0.244 



Table 2f. Detailed Simulation Results (continued) 



Absorber Intercooling: 3''c CO,. 8 
33.5 wt. % MEA,„„ HSS Loading 



^ absorber/-^ base 



Optimum L/G, Cooling at Absoroer Middle 



modeled with a constant temperature (50°C) HEATER, and 
the reflux is fed above segment 1. 

The solvent loop was not closed. It was kept open at the 
absorber inlet. In this way the lean solvent becomes an input 
to the model. Temperature, MEA concentration, and CO2 
loading of the stream are specified. The lean solvent rate is 
calculated with a design specification on the absorber block 
that keeps the removal at a specified value. 

Sensitivity analyses were done on solvent rate, column 
heights, CO2 removal, HSS loading, stripper pressure, and 
absorber temperature. The results of most runs are summa- 
rized in Table 2. 
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in Loading j 



0.15^ 



Absorber (L/G) 



s 2. Optimization of solvent flow rate: 3% mole CO^ 
in flue gas, 85% removal, 0.1 mol HSS/mol 
^AEA^g^, solvent rate normalized to flue-gas 
rate. 



Solvent rate 

In order to optimize the solvent rate, the lean loading was 
varied and the solvent rate was calculated to keep the re- 
moval constant at 85%, with 3% CO. in the flue gas. Figure 
2 demonstrates the optimum. At very knv swlvent rates, ex- 
cessively low lean loading is necessary lo maintain absorber 
performance. At the optimum solvent rate, the specified re- 
moval is achieved with a higher lean loading, which reduces 
the energy required to regenerate the solvent. At very high 
liquid rates, the energy required to heat the solvent to strip- 
per temperature dominates the reboiler duty. 

The heat provided in the reboiler by condensing steam is 
used in the stripper by three energy sinks: the heat of vapor- 
ization of water, the heat of desorption of COj, and the sen- 
sible heat to bring the solvent to reboiler temperature 



Qr = 



+ (F-«H,o)Aiy™p.H.o 

+ Lc^iT,.,^m-T^,) (13) 



§ 1.5 











- 




Operating Line / \ 




















Equilibrium Lin^ 









0.15 0.2 0.25 0.3 0.35 0.4 0.45 
CO^ loading (mol CO^/mol MEA^^^ 

Figure 4. Absorber IVIcCabe-Thiele diagram for opti- 
mum L/G case (L/G^^^ = 0.86): 3% CO2 in 
flue gas, 85% removal, 33.5 wt % IVIEA, 0.1 mol 
HSS/mol MEA^,. 

Figure 3 shows the allocation of the reboiler duty to the three 
terms of Eq. 13. The RateFrac model uses rigorous estimates 
of these thermodynamic quantities. We have used approxi- 
mations to develop Figure 3. Sensible heat was calculated as- 
suming constant L and in the stripper. At the average 
between stripper top and bottom conditions, A//^^^ was as- 
sumed constant. Stripping steam was calculated by the differ- 
ence between the calculated reboiler duty and the first two 
terms. 

Sensible heat is linear, because it is directly proportional to 
L. Stripping steam is nonlinear and tends to an asymptote at 
high L/G. The relative energy requirement for COj desorp- 
tion is constant. 

Optimum Case. Figures 4 and 5 show McCable-Thiele di- 
agrams for absorber and stripper at the optimum solvent rate. 



5' 











Q 1 

strippin 


g Steam'" 


C02 












Q 


Q 


sensible l^t C02 




1 V i \ 



Figure 3. Allocation of the reboiler duty to CO2 desorp- 
tion, sensible heat, and stripping steam gen- 
eration. 















Eauilibrium Line • / 






/ / 




y '^jr OpevaimQ Line 


- ,y 









0.15 0.2 0.25 0.3 0.35 0.4 0.45 
CO^ loading (mol CO^/moi MEA^^^) 

Figure 5. Stripper IVIcCabe-Thiele diagram for optimum 
L/G case (L/G^^^ = 0.86): 3% COj in flue gas, 
85% removal, 33.5 wt % MEA, 0.1 mol 
HSS/mol MEA,„,. 
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/ 






Equilibrium Line J 












/' 

/ y 

^*'*Operat 












ng Line 













0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 
CO^ loading (mol CO^/mol MEA^^^^) 

Figure 6. Stripper McCabe-Thlele diagram for low L/G 
case (L/Q^ = 0.63): 3% CO2 In flue gas, 85% 
removal, 33.5 wt. % MEA, 0.1 mol HSS/mol 
MEA^,. 




0.26 0.28 0.3 0.32 0.34 0.36 0.38 0.4 
CO^ loading (mol CO^/mol IVIEA^q^ 

Figure 7. Stripper McCabe-Thlele diagram for high 
L/G case (L/G„3^ = 1.77): 3% CO2 In flue 
gas, 85% removal, 33.5 wt. % MEA, 0.1 mol 
HSS/mol MEA^,. 



The operating lines are obtained from the simulation. The 
circles drawn on them represent each RateFrac segment, with 
the same height of packing. For this sensitivity analysis, 20 
segments were used lo represent the absorber packing and 30 
to represent the shipper packing. The equihbrium lines are 
obtained from llash calculations at the temperature and com- 
position of three selected segments. 

At the optimum solvent rate a rich pinch is obsei-ved in the 
stripper. Nevertheless, the separation is well distributed 
among the 30 segments. The rich pinch is a consequence of 
the nonlinear operating line. Water condenses from the va- 
por going up the column, thus curving the operating line to- 
ward the equilibrium line. The fact that the top of the strip- 
per is approximately 10°C colder than the bottom reduces the 
nonlinearity of the equilibrium line, contributing to the for- 
mation of the rich pinch. A considerable separation is 
achieved in the reboiler stage, because it was modeled as an 
equilibrium segment. 

At the rich end of the stripper there is a discontinuity in 
the loading at the feed stage. This is a consequence of flash- 
ing the hot feed to the lower equilibriimi temperature at the 
top of the stripper. 

The absorber McCabe-Thiele diagram shows that no equi- 
librium pinches are present, and a wide driving force is well 
distributed through all the segments, with a slightly higher 
concentration of segments at low loading, indicating that the 
optimum case is closer to a lean pinch rather than a rich 
pinch. 

Low L/G. Figure 6 shows the stripper McCabe-Thiele 
diagram with a low solvent rate (L/G^^^^ = 0.63). The lean 
loading is low (0.11), in order to maintain the removal at 85% 
with the lower solvent rate. A tight pinch is observed in the 
middle of the stripper. The nonlinearity of the equilibrium 
relationship requires that the equilibrium be approached at 
an intermediate loading. Figure 3 shows that a large amount 
of stripping steam is required to overcome the middle pinch. 

In the absorber there were no evident pinches, even though 
a rich pinch would be expected at low L/G. 



High L/G. Figure 7 shows the McCabe-Thiele diagrams 
for the stripper at a high solvent rate. The stripper rich pinch 
is more definite, with more than half of the column working 
at the rich end, pinched conditions. From Figure 3, it can be 
seen that the stripping steam contribution to the heat re- 
quirement does not change significantly at high L/G. The 
sensible heat determines an increase in reboiler duly. The 
feed loading effect is reversed. Flashing is reduced, because 
of the lower rich loading. In this case, the loading increases 
at the top of the stripper, and it keeps increasing in the top 
half of the column. This is associated with CO 2 absorption in 
the top half of the column. The high liquid rate generates 
high liquid heat capacity, the vapor cannot heat up the liquid 
feed fast enough, and it gets cooled by it. Water condensa- 
tion occurs, causing the CO 2 to become more concentrated 
in the vapor phase, to the point where absorption occurs. 

There is no apparent equilibrium pinch in the absorber. 
However, half of the column is operating between 0.27 and 
0.3 loading, suggesting a rate-based "pinch" at the lean end. 

Effects of absorber height, CO2 concentration, and CO 2 
removal 

The absorber height was varied about the value of the base 
case, with 3% C02/85% removal, 3% CO2/90% removal, and 
10% CO2/90% removal (Figure 8). The absorber diameter 
and the stripper size were kept constant at the base -case 
value. The solvent rate was optimized for each case. The re- 
boiler duty decreases with packing height, tending to an 
asymptote. The minimum duty was calculated increasing the 
height of the absorber at the point where the reboiler duty 
stopped changing. For the three curves, the values of the 
mimimum reboiler duties (relative to the reference arbitrary 
value) are 0.9425 (3% CO2, 85% removal), 0.945 (3% CO2, 
90% removal), and 0.891 (10% CO2, 90% removal). The min- 
imum reboiler duty is practically independent of CO2 re- 
moval. 
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Figure 8. Effect of absorber height on reboiler duty per 
mole of CO2 removed, 33.5 wt. % lUlEA, 0.1 mol 
IHSS/mol IVIEAt^, optimum solvent rate. 

With a finite amount of packing, the reboiler duty in- 
creases with CO2 removal, because the optimum solvent rate 
increases as the removal increases. This increases the sensi- 
ble heat requirement. 

At constant removal, the normalized reboiler duty in- 
creases as the CO2 concentration in the flue gas decreases. 
Even though the absolute reboiler duty is higher for the 10% 
CO2 case than the 3% COj case, the reboiler duty normal- 
ized to the total number of moles recovered is lower. A sim- 
ple explanation can be foimd if the minimimi thermodynamic 
work required for the process is calculated at isothermal con- 
ditions (Eq. 14) 

Pen out 

^-^ = AG,. ,, = /;rin ^ ^' (14) 

I'or the case of 3 mol % of CO 2 in the flue gas, Pcoj.in = 
0.03 atm; for the 10% case, P^o^ = 0.1 atm. Assuming an out- 
let pressure of 2 atm, the minimum work required to com- 
press the CO 2 from the absorber inlet condition to the strip- 
per outlet condition is ln(2/0.03)/ln(2/0.1) = 1.4 times higher 
in the low COj case. 

With 3% CO2 and 85% removal at ^absorberAbase case = 1-6, 
the reboiler duty is practically at its asymptotic value. A rich 
pinch would be expected in the absorber. Figure 9 shows the 
absorber McCabe-Thiele diagram for this case. The ab- 
sorber does not appear to have an equilibrium pinch. How- 
ever, the distribution of packing segments, represented by the 
circles on the operating line, shows that several segments 
achieve little change at the rich end. This apparent rate -based 
pinch is not determined by a low driving force. At the rich 
conditions there is a reduced concentration of free amine to 
enhance the reaction rates, therefore, the very low reaction 
rates reduce the rate of mass transfer. 

Stripper height 

The stripper height was varied to values less than the base 
case, for 3% COj and 85% removal (Table 2c). In this analy- 
sis the liquid-phase mass-transfer performance in the stripper 




0 Lrr-l^a , \ i ! 
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CO^ loading (mol COj^/mol MEA^^^) 



Figure 9. McCabe Thiele diagram for the absorber: 

^absorber/^base case = 1 -6, 33.5 Wt % IWEA, 0.1 mOl 

HSS/mol MEA,o,, optimized solvent rate 
(L/G„,,, = 0.82), 3 mol. % COa in flue gas, 
85% removal. 



is overpredicted by the model, due to the use of instanta- 
neous reactions. The stripper performance will also be over- 
predicted because the reboiler is modeled as an equilibrium 
segment. 

The importance of the error introduced by using instanta- 
neous reactions was estimated calculating gas and liquid 
mass-transfer resistances. The gas film resistance ranges from 
50% to 70% in the base-case stripper. Since the reaction-rate 
resistance represents 20% to 70% of the liquid resistance 
(Freguia, 2002), 10% to 35% of the total mass-transfer resis- 
tance is neglected in the model. 

The stripper height affects the heat requirement less dra- 
matically than the absorber height. A column half as high 
could perform almost as well, with an increase in reboiler 
duty less than 4%. 

This conclusion may be valid only for 3% COj. Figure 10 
shows the McCabe-Thiele diagram for the stripper at 10% 
CO 2 at the optimum solvent rate. The stripper does not have 
an equilibrium pinch, as with 3% CO2. This is because there 
is 3.3 times more CO2 to transfer through the same gas-liquid 
contact area. Therefore, the stripper height should have a 
larger effect on reboiler duty with 10% COj. 

Heat stable salts 

The base case was calculated with HSS loading of 0.1 to 
represent probable effects of MEA degradation. The partial 
neutralization of MEA with a strong acid can increase the 
equilibrium partial pressure of COj at a given temperature 
and loading, because strong acids shift the equilibrium to- 
ward the acid species CO 2- Figure 11 shows the McCabe- 
Thiele diagram for the stripper with no HSS and with 0.1 mol 
HSS/mol MEA^t, at fixed solvent rate and lean loading. A 
reduction of 40% in the reboiler duty is obtained upon adding 
10% HSS. With HSS loading = 0.1, the operating line follows 
the equilibrium line, keeping the driving force almost con- 
stant, typical of an optimized process. Without HSS, the re- 
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Figure 10. IVIcCabe-Tliiele diagram for the stripper: 

^absorber/^base case = 1 . 33-5 % MEA, 0.1 mol 

HSS/mol MEA,„t, optimized solvent rate 
(VGrnass = 2.68), 10 mol % COa in flue gas, 
90% removal. 

boiler duty is 40% higher, because the solvent rate is not 
optimized. The addition of HSS changes the optimum solvent 
rate and the optimum lean loading. 

The positive effects of HSS have been obsei-ved in other 
absoiption/stripping systems with nonlinear equilibrium rela- 
tionships (Kohl and Nielsen, 1997; Carey, 1990). 

When the solvent rate is optimized, the effects of HSS are 
more limited. In Figure 12 the effect of HSS on reboiler duty 
is shown for three amine concentrations, with solvent rate 
optimized for each case. The effects on optimum lean load- 
ing and L/G are shown in Table 2d. The neutralization of 
10% of the MEA reduces the reboiler duty by approximately 
2%, whereas at nonoptimized conditions the reduction was 
40%. 




.2 0.24 0.28 0.32 
CO^ loading (mol CO^/mol MEA^^j.^) 

Figure 11. Comparison of base case Aspen Plus runs 
with 0 and 0.1 mol HSS/mol MEAt.^, heat- 
stable salts added on top of 30 wt. % active 
MEA: 3.13% COj in flue gas, L/G„,,, = 0.78. 




0.05 0.1 0.15 
HSS loading ([HSS]/[MEA]^^^) 

Figure 12. Effect of MEA concentration and of HSS 
loading on the reboiler duty, 3 mol % COj, 
85% removal, optimum solvent rate. 



Reboiler duty is reduced by increasing amine concentra- 
tion, following Eq. 8. However, amine concentrations greater 
than 36% were not studied, because the diffusivity of CO 2 
becomes so low that the absorption rates start to decrease, 
according to Eq. 8, and because MEA solvents become more 
corrosive when concentrated. 

With 30 wt % MEA the minimum reboiler duty occurs at 
HSS loading of 0.2, indicating that HSS can slightly improve 
the performance of the process. Some HSS serves to lin- 
earize the equilibrium relationship and reduce steam require- 
ment. However, at greater HSS, a high solvent rate is neces- 
sary, which increases reboiler duty because of sensible heat. 

Effect of stripper pressure 

Historically, stripper pressure has been designed for 1 to 2 
atm to vent CO 2 and to minimize amine degradation, which 
can occur at a higher pressure/temperature (Rochelle et al., 
2001). Because the ratio Pco^/Pyi^o increases significantly 
with temperature, the reboiler heat duty will decrease with 
increased stripper pressure. However, in a thermally inte- 
grated power plant, it may be attractive to operate at lower 
pressure, because of the reduced cost of using lower pressure 
steam to drive the reboiler. 

The effect of stripper pressure was analyzed at constant 
removal (85%) and CO2 concentration in the flue gas (3%). 
The solvent rate was optimized, giving at P = 0.5 atm, 
V^mass = 1-93; and at P = 3.7 atm, L/G^ = 0.66. Figure 
13 gives the energy requirement for stripping as a function of 
the stripper pressure. Figure 14 gives the effect of pressure 
on the McCabe-Thiele diagram. The reboiler duty increases 
significantly at lower pressure, but only decreases slightly 
above 2 atm. The asymptote represents a minimum reboiler 
dut)', when there is no stripping steam generated and the va- 
por is practically pure COj. The lost work of the reboiler 
steam was calculated with a Carnot efficiency, using a sink 
temperature of 40°C and a loss of ITC for heat transfer across 
the reboiler. The work value of the reboiler steam has a mini- 
mum at about 1.5 atm, and only increases by approximately 
30% at 0.5 atm. Therefore, the stripper can operate effi- 
ciently at any convenient level of steam pressure available 
from the power plant. The calculated lost work does not in- 
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Figure 13. Effect of stripper pressure on the reboller 
duty and lost worl<, for a 3 mol % COg flue 
gas and 85% removal, 33.5 wt. % MEA, 0.1 
mol HSS/mol IVIEA,^,, L/G optimized, 20°F 
approach in cross-exchanger. 



elude the additional compression cost for the COj, but that 
is small relative to the value of the reboiler steam. Table 2e 
presents more detailed results of this sensitivity analysis. 

For a solvent with the ratio - ^H^i^^ qq^AH^^^^^^q equal 
to 1, the stripper pressure would not affect the reboiler duty. 
The stripper could be run at low pressure with lower-lost 
work. The compression cost would be part of the pressure 
optimization. Another advantage would be a lower amine 
degradation rate, due to lower operating temperature. Exten- 
sive work is underway to study solvents with a lower heat of 
absorption than MEA, such as piperazine -promoted KjCOj 
((\illinanc, 2002). 

Effect of intercooUng 

Figure 15 shows that the temperature in the middle of the 
absorber can be 15-35°C hotter than the temperature of the 
flue gas or solvent at either end. This phenomenon is caused 
by the heat of absorption of CO 2 and is more significant with 
10% CO2 than with 3% CO2. The absorption rate of CO2 
can increase with temperature to the extent that the com- 
bined rate constant 



^/c2[MEA]Z)co, 

increases with temperature. The rate of absorption will de- 
crease to the extent that it is determined by the equilibriimi 
driving force, Pqo^ — Poa^' which decreases with tempera- 
ture. 

Figure 15 shows that, for a 10% CO; absorber, the inter- 
cooling reduced the hquid temperature to 42°C in the middle 
of the absorber. Intercooling has practically no effect on the 
temperature in the rich end of the column. An amount of 
heat equal to one-third of the reboiler duty was removed from 
the liquid at the middle point of the absorber. The solvent 
rate was optimized and the removal was kept constant at 90%. 
However the reboiler duty was reduced only by 3.8%. Figure 
16 compares the absorber McCabe-Thiele diagrams for the 
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Figure 14. McCabe Thiele diagrams for strippers oper- 
ating at 3.7 and 0.5 bar, 3 mol % COj, 85% 
removal, 33.5 wt. % MEA, 0.1 mol HSS/mol 
MEAtot, optimized L/G. 

cooled and noncooled cases. The equilibrium line is lower for 
the cooled case, increasing the driving forces throughout the 
column, but only slightly improving the absorber perfor- 
mance. 

Intercooling with 3% CO2 in the flue gas will be less effec- 
tive, because the temperature bulge is less significant. 

Conclusions 

With optimized solvent rate at the base-case conditions, 
the stripper approaches equilibrium at the rich end and has a 
well-distributed driving force throughout the rest of the con- 




Figure 15. Comparison of liquid and gas temperature 
profiles for 10% CO2 90% removal inter- 
cooled absorber, 10% CO2-90% removal 
nonlntercooled absorber, and 3% C02-85% 
removal nonlntercooled absorber: Inter- 
cooler at absorber middle point, cooling duty 
1/3 reboiler duty, 0.1 mol HSS/mol MEA, 
optimized L/G. 
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-no intercooling, 10% CO ; 
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Figure 16. Comparison of absorber IVIcCabe-Thiele di- 
agrams for jntercooled and nonintercooled 
absorber, intercooler at half column, inter- 
cooler duty 1/3 reboiler duty, 10% COg, 90% 
removal, 0.1 mol HSS/mol IVIEA, optimized 
L/G. 

tactor. The rich pinch results from the decreased vapor rate 
at the top of the column. The absorber does not have a close 
approach to equilibrium, with high driving forces throughout 
the column. 

The solvent cireulalkni rale alTeets the steam requirements 
heavily, and can be optimized I'or each set of design condi- 
tions. A nonoptimized process can consume 50% more en- 
ergy if the solvent rate is 30% smaller than the optimum, or 
25% more if it is twice the optimum. The loss in performance 
is more important if the solvent rate is low, therefore, it is 
preferable to operate the process at a solvent rate slightly 
higher than the optimum. All of the following conclusions 
assume an optimized solvent rate. 

Steam consumption is reduced by a greater height of pack- 
ing in the absorber. I?y providing more area for mass trans- 
fer, it is possible to enhance the process performance. I'hc 
effect is important with less packing, whereas with a large 
amount of packing the effect disappears and the reboiler duty 
reaches an asymptote, only 10% less than the base-case de- 
sign. Even with a considerable height of packing, the ab- 
sorber does not approach equilibrium. In reality, at the rich 
end many stages achieve very little removal. The "pinch" is 
determined by slow kinetics, rather than low driving forces. 

Steam consumption is only reduced slightly by increased 
packing height in the stripper, at the base -case conditions of 
this study. The stripper height required is much less than that 
of the absorber, because, due to higher temperature, the re- 
action rate is not limiting. The amount of packing in the 
stripper should be scaled according to the amount of COj 
removed. 

An increase in COj content in the flue gas and the reduc- 
tion of COj removal are accompanied by a small reduction in 
the specific steam consumption (kJ/mol COj). The effects 
tend to disappear if the columns are oversized. 

The addition of a strong acid to the solvent forms HSS, 
which can slightly reduce the energy requirement, by making 



the equilibrium relationship more linear. At a constant sol- 
vent rate, the addition of HSS can change steam consump- 
tion by as much as 40%. 

The stripper pressure affects the ratio of C02:H20 in the 
stripper vapors, affecting the reboiler duty and the lost work. 
An optimum pressure is found at approximately 1.5-2 atm. 
The effect is not very significant, therefore, it is possible to 

Intercooling increases absorption rates in the absorber, by 
increasing the driving force for absorption. The effect is mod- 
erate, since a heat stream equal to one-third of the reboiler 
duty reduces the latter by only 4%, when extracted from the 
center of the absorber. 
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Notation 

Ocoj,; = activity of CO2 at the gas-liquid interface 
Ocoji = activity of CO2 in equilibrium with the interface coi 
position 

^tot = tot^l concentration of liquid phase, M 
D = column diameter, m 
Deo =diffusivity of CO, in solution, m^s 
I»co„wate?=diffusivity of CO2 in water, m^s 
HSS loading =heat stable salts loading 

G =gas flow rate in absorber, kmol/s 
^^abs,co2 =heat of absorption, kJ/mol 
A//^3pj^^Q =heat of vaporization of H2O, kJ/mol 
Hco-i = Henry's constant of CO2, atm/M 
[ /]] = concentration of species ; at the interface, kmol/m' 
/f , = rate constant of second-order reactions, m ^kmol s 

equilibrium constant 
L = liquid flow, kmol/s 
ATco, =flux of COo, kmol/m^s 
Hco \ = moles of CO, removed from flue gas, kmol/s 
P = total pressure, atm 
(?r = reboiler duty, MJ 
O^^j = relative reboiler duty, dimensionless 
i?= reaction rate, M/s 
R = universal gas constant, 8.314 kJ/kmol K 
r = temperature, K 
K= vapor flow rate in stripper, kmol/fe 
Xj = liquid mole fraction of species ; 
y: =gas mole fraction of species j 
Z = height of packing in a colimm, m 



Greek letters 



a =C02 loading 

y, = activity coefficient of species j with reference si 
pure component 
=viscosity of solution 
=viscosity of water 
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